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INTRODUCTION 

An experimental study of mixing downstream of axial and swirling coaxial jets 
Is being conducted to, obtain data for the evaluation and Improvement of turbulent 
transport models currently employed In a variety of computational procedures used 
throughout the propulsion community. The nonswirling coaxial jet study was 
completed under Phase I and is reported In Reference I. The swirling coaxial jet 
study was conducted under Phase II of the contract and Is reported In Reference 2. 

A TEACH code was acquired, checked out for several test cases, and Is reported In 
Reference 3. A study to measure length scales and to obtain a limited number of 
measurements with a blunt trailing edge Inlet Is being conducted under Phase I I I of 
the contract and HOST sponsorship. Results from Phases I and II have also been 
presented in References 4, 5 and 6. 


PHASE III OBJECTIVES 

The Phase III effort was directed toward (I) the aqulsltlon of length scale and 
dissipation rate data that will provide more accurate inlet boundary conditions for 
the computational procedures and a data base to evaluate the turbulent transport 
models In the near jet region where recirculation does not occur and (2) the 
acquisition of mass and momentum turbulent transport data with a blunt inner- jet 
inlet configuration rather than the tapered inner- jet Inlet employed In Phases I and 
I I. 


Mass and momentum turbulent transport data was obtained downstream of the blunt 
Inner-jet inlet configuration, using the laser vel oc Imeter/ I aser induced f I uoresence 
measurement techniques employed In the previous experiments with swirling and 
nonswirling flow conditions. The velocities, concentrations, mass turbulent 
transport and momentum turbulent transport results obtained downstream of the blunt 
inner- jet inlet configuration for the nonswirling flow condition were not 
significantly different from those obtained with the tapered Inner-Jet Inlet. The 
results will be presented and compared with the Reference 1 results in the 
forthcoming Phase III Interim summary report. 

The measurement technique generally used to obtain approximate integral length 
and micro scales of turbulence and dissipation rates was recently computerized at 
UTRC. This computerized data aqulsltlon and reduction procedure was then used to 
obtain length scale and dissipation rate data for three test configurations. The 
data and results obtained with the tapered Inner jet Inlet configuration and the 
nonswirling flow condition will be discussed In the following sections. 


^Performed under Contract NAS3-22771 . 
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LENGTH SCALE AND DISSIPATION RATE MEASUREMENTS 


Data. Analysis Formulation 


With Taylor’s hypothesis, the temporal variation of a one dimensional velocity 
fluctuation Is related to a spatial variation by the convective velocity, U. This 
Is also referred to as the frozen-turbulence approximation (Ref. 7). The auto- 
correlation is defined 

R(r ) = u ( t ) u(t+T )/u' 2 

An Integral scale and an approximate micro scale of turbulence and an eddy 
dissipation rate are related to the autocorrelation. The integral scale of 
turbulence Is: 


LI = U J o °°R (r ) dr 


A second Integral scale of turbulence Is also obtained directly from spectrum 
(Ref. 7) 


L2 = ( U/4 ) I im n _, 0 [E(n)/u’ 2 ] 

The microscale length of turbulence, A, Is defined as the Intercept of the parabola 
that matches R(r) at r = 0, I.e., 

R( t ) = 1 - (tU/A) 2 

The dissipation rate, £ , Is related to the microscale of turbulence, the rms 
velocity and the fluid kinematic viscosity, : 

e = 30 v u’ 2 /A 2 

Although the autocorrelation/microscale relations shown above will lead to a 
dissipation rate, the determination of the curvature of R(r) at t = 0 is often 
difficult and sometimes ambiguous. Another relationship for the dissipation rate is 
obtained directly from the one-dimensional velocity turbulent energy spectra 

1 /A 2 = ( 27 T)/(UV 2 ) r n 2 E ( n ) dn 

,J o 

where E(n) is the turbulent kinetic energy per frequency cycle of one velocity 
component at frequency n. Thus the approximate dissipation rate is more easily 
obtained with an independent spectrum measurement. 

An autocorrelation for the spectral data is also constructed from the Fourier 
transform of the energy spectral distribution using a computer routine which Is part 
of the data reduction program. The Fourier transform which produces the spatial 
correlation is 

R (T ) = rv™ E(u) d(j 

• ; o 


The information flow chart for the computerized length scale measurement technique 
Is shown in Figure 1. 
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Je-5 1. Cond 1 tions anxL Apparatus 

The length scale and dissipation rate measurements were conducted In the UTRC coaxial 
flow facility (Figure 2) using water as the working fluid. The flow conditions and 
geometries were Identical to those employed in Reference 1, i.e. U (inner) = 0.52 
m/s and U (annular) = 1.66 m/sec and with a Reynolds number of 35,000 based on the 
duct diameter and the average velocity. 

A sketch of the Inlet and measurement locations is shown In Figure 3. The 
swirler was not Installed for the measurements discussed herein. Mean, fluctuating 
and spectral measurements of the axial velocity were obtained In the regions without 
recirculation at z = -41, 5, 51 and 102 mm from the inlet plane. Mass and momentum 
turbulent transport data was obtained at the latter three locations during Phase I 
(Reference 1). Standard hot film constant resistance measurement procedures were 
used to aquire the data. The water was filtered to prevent particles in the water 
from attaching to the probe and altering the probe calibration curve. 


Results 

The mean and fluctuating axial velocity measurements obtaining at the four 
axial locations are presented In Figures 4 and 5 and compared with results from 
Reference 1 obtained with a laser ve loci meter. The agreement between the two sets 
of measurements Is excellent. The conclusion from this comparison is that the hot 
film data aquisition and reduction procedure for obtaining mean and fluctuating 
velocities reproduces the previous results with acceptable accuracy. 

The dissipation rate distributions for the four axial locations are presented 
in Figure 6. Note the dissipation rates at the center of the inner jet and the 
annular jets are approximately the same at z = -41 and 5 mm from the inlet plane. 

At z = 5.1 mm, the dissipation rate across the inlet varies by 2 1/2 orders of 
magnitude. This radial variation decreases to less than one order of magnitude (by 
Increasing the lower dissipation rates and decreasing the higher dissipation rates) 
at z = 102 mm from the inlet plane. As expected the locally high and low 
dissipation rates are associated with high and zero axial velocity shear rates, 
respectively. 

The integral length and micro scales of turbulence are presented In Figure 7. 
The integral scales LI and L2, defined on the figure, were both determined as part 
of the computerized data reduction procedure. Note that the integral scales 
determined by both methods are In generally good agreement. The largest difference 
between the two scales is about 20 percent and occured at z = 102 mm. This is a 
location with long wave length eddies occurlng over a moderate frequency range. The 
differences between the two Integral length scales are generally less than ten 
percent with the Integral scale determined from the integration of the auto 
correlation always being greater. The integral scales of turbulence in the inner 
jet inlet were 9 to 12 mm, approximately half the tube diameter and compatible with 
length scales reported In Reference 7 for previous work by Laufer. 

The micro seal es of turbulence ,A, were 15 to 40 percent of the integral length 
scales. The ratio of micro scale to integral scale inside the inlet duct (z = -41 
mm) was somewhat greater than measured by Laufer for a two-dimensional duct, i.e. : 
15 percent of the UTRC circular duct diameter compared to 10 percent of the Laufer 
2-d slot height. Thus, the present results in the duct appear to be reasonably 
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compatable with previous studies. The micro scale of turbulence does not vary much 
( 1 < x < 3.8mm) compared to the Integral scale of turbulence. For the two 
measurement locations downstream of the Inlet, the ratio of the maximum to minimum 
microscales (2<X<4mm) Is only half the ratio of the Integral length scales 
(4.2<L <1 9) . 
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FIG. 1 


INFORMATION FLOW CHART FOR COMPUTERIZED 
LENGTH SCALE MEASUREMENT TECHNIQUE 
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MEAN AXIAL VELOCITY, U - m/sec 


MEAN AXIAL VELOCITY PROFILES FOR NONSWIRLING FLOW 
WITH TAPERED INNER JET INLET 
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FLUCTUATING AXIAL VELOCITY, u '-m/sec 


FLUCTUATING AXIAL VELOCITY PROFILES FOR NONSWIRLING FLOW 
WITH TAPERED INNER JET INLET 
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